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We study the problem of computing the probability for the time-of-arrival of a quantum particle at 
a given spatial position. We consider a solution to this problem based on the spectral decomposition 
of the particle's (Heisenberg) state into the eigenstates of a suitable operator, which we denote as 
the "time-of-arrival" operator. We discuss the general properties of this operator. We construct the 
operator explicitly in the simple case of a free nonrelativistic particle, and compare the probabilities 
it yields with the ones estimated indirectly in terms of the flux of the Schrodinger current. We derive 

os : 

a well defined uncertainty relation between time-of-arrival and energy; this result shows that the well 
known arguments against the existence of such a relation can be circumvented. Finally, we define 
a "time-representation" of the quantum mechanics of a free particle, in which the time-of-arrival 
is diagonal. Our results suggest that, contrary to what is commonly assumed, quantum mechan- 
ics exhibits a hidden equivalence between independent (time) and dependent (position) variables, 
analogous to the one revealed by the parametrized formalism in classical mechanics. 



I. INTRODUCTION 

Consider the following experimental arrangement. A particle moves in one dimension, along the x axis. A detector 
is placed in the position x — X. Let T be the time at which the particle is detected, which we denote as the 
"time-of-arrival" of the particle at X. Can we predict T from the knowledge of the initial state of the particle? 

In classical mechanics, the answer is simple. Let x(t] xq,po) be the general solution of the equations of motion 
corresponding to initial position and momentum xo and po at t = 0. We obtain the time-of-arrival T as follows. We 
invert the function x — x(t;xo,po) with respect to t, obtaining the function t(x;xo,po). The time of arrival T at X 
of a particle with initial data xq and po is then 



a: 



T = t(X;x ,p ). (1) 



2 . 

Two remarks are in order. First, if t(x;xo,Po) is multivalued, we are only interested in its lowest value, since the 
• • i particle is detected the first time it gets to X. Second, for certain values of Xq and po, it may happen that X is 
outside the range of the function x(t; Xq,Pq). This indicates that the detector in X will never detect a particle with 
that initial state. The time-of- arrival is a physical variable that -in a sense- can take two kinds of values: either a 
real number T G H, or the value: "T = newer" . Notice that in the latter case, the quantity T formally computed 
from ([j]) turns out to be complex. Thus, a complex T from (Q) (for given X;xq,Po) means that the particle with 
initial data xo,po is never detected at X. 

In quantum mechanics, the problem is surprisingly harder. In this case the time of arrival can be determined only 
probabilistically. Let tt(T) be the probability density that the particle is detected at time T. Namely, let 

T 2 

n{T)dT (2) 

Ti 

be the probability that the particle is detected between the time T\ and the time Ti. How can we compute 7r(T) from 
the quantum state, e.g. from the particle's wave function ip(x) at t = 0? 
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To the best of our knowledge, this question has not received a complete treatment in the standard literature on 
quantum mechanics. The problem of computing the time of detection of a particle is usually treated in very indirect 
manners. For instance, the probability (in time) of detecting a decay product - say a particle escaping the potential 
of a nucleus- can be obtained from the time evolution of the probability that the particle is still within the confining 
potential. Alternatively, one can treat the detector that measures the time-of-arrival quantum-mechanically, and 
compute the probabilities for the positions of the detector pointer at a later time; in this way one can trade a 
meaurement of the time-of-arrival for a measurement of position at a late fixed time. In the fifties, Wigner considered 
the problem of relating the energy derivative of the wave function's phase shift to the scattering delay of a particle 
0. This approach, later developed by Smith [g] and others (see for instance Gurjoy and Coon ||) gives the average 
delay, but fails to provide the full probability distribution of the time of arrival. (Smith's paper begins with: "It is 
surprising that the current apparatus of quantum mechanics does not include a simple representation for so eminently 
observable a quantity as the lifetime of metastable entities.") In the seventies, Piron discussed the problem in a 
conference proceeding Q, sketching ideas related to the ones developed here. Ideas related to the ones presented 
here were explored in ref. g, but in this case too only the average time-of-arrival was obtained, and not its full 
probability distribution. Kumar || studied the quantum first passage problem in a path integral approach, but did 
not obtain a positive probability density. The problem has been studied in the framework of Hartle's generalized 
quantum mechanics |7]] by using sum over histories methods. Various attempts in this direction and discussions of 
difficulties can be found in Ref. || . See also the recent paper ref. |^] for a discussion of the problem and for references; 
in particular, ref. Q discusses the difficulties one has to face in trying to compute sequences of times-of-arrival - 
an important problem which, however, we do not address here. As stressed by Hartle in |9), generalized quantum 
mechanics generalizes "usual quantum mechanics" ; here, on the other hand, we are interested in the question whether 
7r(T) can be computed within the mathematical framework of conventional Hamiltonian quantum mechanics. 

We see two reasons of interest for discussing the problem of computing time-of-arrival in quantum mechanics. First, 
it is a well posed problem in simple quantum theory, and there must be a solution. Echoing Smith 0, we do not 
expect that quantum mechanics could fail to predict a probability distribution that can be experimentally measured 
by simply placing a detector at a fixed position and noting the time at which it "clicks" . The problem is not just 
academic: it is related to the problem of computing the full probability distribution (as opposed to the expectation 
value) for the tunnelling time through a potential barrier. This problem has relevance, for instance, in computing 
rates of chemical reactions (see, for example, Kumar ||). Second, the problem bears directly on the interpretation 
of quantum theories without Newtonian time [ fbf|]l5f| and thus on quantum gravity; we shall briefly comment on this 
issue in closing. 

This paper is the first of a sequence of two. Here we develop a general theory for the time-of-arrival operator, and 
study the free nonrelativistic particle case in detail. In a companion paper pO| j, we investigate a technique for the 
explicit construction of the time-of-arrival operator in more general cases, we extent our formalism to parametrized 
systems, and we study some less trivial models: a particle in an exponential potential and a cosmological model. 

In the next section we give a general argument, based on the superposition principle, for the existence of an 
operator T (the time-of-arrival operator) such that tt(T) can be obtained from the spectral decomposition of ip(x) 
in eigenstates of T, in the usual manner in which probability distribution are obtained in quantum theory. T has 
peculiar properties that distinguish it from conventional quantum observables. We give a general argument based on 
the correspondence principle indicating that T can be expressed in terms of position and momentum operators by the 
inverse of the classical equations of motion, eq. (Q). This does not suffice in fixing the operator, since factor-ordering 
ambiguities can be serious. The problem of the actual construction of the operator T in more general systems will 
be addressed in []l0| . In Section 3 we study an explicit form of the operator in the case of a free nonrelativistic 
particle. We diagonalize the operator, providing a general expression for the time-of-arrival probability density 7r(T). 
In particular, we calculate 7r(T) explicitly for a Gaussian wave packet. In section 4 we discuss some consequences of 
our construction. We notice that the existence of the operator implies that the quantum mechanics of a free particle 
can be expressed in a "time-representation" basis. We derive time-energy uncertainty relations. We conclude in 
section 5 with a general comment on the equivalence between time and position variables suggested by our results. 

In the Appendix, we study whether the probability distribution we computed is reasonable, by comparing it with 
the one estimated indirectly using the Schrodingcr current. We find that the two agree within second order in the 
deBroglie wavelength of the particle. The probability computed from the Schrodinger current cannot be physically 
correct to all orders because it is not positive definite; whether or not the probability distribution computed with T 
is physically correct to all orders is a question that can, perhaps, be decided experimentally. 



II. TIME-OF-ARRIVAL: GENERAL THEORY 



2 



A. The incomplete spectral family P(T) 



Consider the quantum analog of the experimental situation sketched at the beginning of the paper: a particle is 
in an initial state ip at t = 0, and a particle detector is placed at x = X. Let T be the time at which the particle is 
detected. Let TV^(T)dT be the probability that the particle is detected between times T and T + dT. Let ip and <j> 
be two quantum states such that both n^(T) and 7i>(T) have support in the interval I = (T, T + AT). Consider the 
state formed as the linear combination aip + b<p (where a and b are any two complex numbers with \a\ 2 + \b\ 2 = 1). 

According to the superposition principle, if a measurable quantity has a definite value when the system is in the 
state ip and value \^ when the system is in the state <p, then a measurement of such a quantity in the state aip + b<p 
will yield either A^, or A^ (with respective probabilities \a\ 2 and \b\ 2 ) [O. If we assume the general validity of the 
superposition principle, we must then expect that the probability distribution TT a ip+b4>(T) will have support in the 
interval (T, T + AT) as well. Therefore, the states ip such that tv^(T) has support on a given interval / = (T, T + AT) 
form a linear subspace of the state space. We can therefore define a projection operator Pj as the projector on such 
a subspace. 

The superposition principle could fail for the time-of-arrival. However, we would be surprised if it did. Notice 
that the question can (probably easily) be decided experimentally. Perhaps an experiment testing the validity of the 
superposition principle in this contest could have some interest. Here, we assume that the principle holds, and thus 
the projectors Pt,t+at are well defined. 

By their very definition, the projectors satisfy P7P7 = Pj if the interval J is contained in the interval /, and 
PjPj = if the two intervals are disjoint. The operators Pj can therefore be written in terms of a family of (spectral) 
projectors P(T) as 

pT+AT 

Pt,t+at = J t P{T')dT'. (3) 

The (spectral) family P(T) contains all the information needed to compute 7r(T). Indeed, from the definition given, 
and using again the superposition principle, we have easily 

^(T) = (i,\P{T)\i,). (4) 

Thus the probability distribution n(T) can be obtained in terms of the spectral family P(T), in the same way in which 
all probability distributions are obtained in quantum mechanics. Indeed, recall that if A is the self-adjoint quantum 
operator corresponding to the observable quantity A, then the probability distribution tt^(A) of measuring the value 
A on the state ip is ( ip \ P(A) \ip), where P(A) is the spectral family associated to A, namely 

A= f A P{A) dA. (5) 



B. The operators T and P 

One may be immediately tempted to define a "time-of-arrival operator" in analogy with (S) as 

r+co 

f=\ TP(T)dT (6) 



so that an eigenstate of this operator with eigenvalue T would be a (generalized) state detected precisely at time 
T. However, there is an important difference from usual self-adjoint quantum mechanical observables that must be 
addressed before doing so. If P(A) is the spectral family of a self-adjoint operator A, then 

P(A) dA = 1, (7) 

-oo 

where 1 is the identity operator. On the other hand, define P by 

>c 

P(T) dT; (8) 
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there is no reason for P to be the identity. If it is not, we say that the spectral family P(T) is "incomplete". 

Incompleteness occurs because it is not true that any state is certainly detected at some time. Most likely, there 
are states that are never detected, given that such states exist in the classical theory as well. Thus, P projects on 
the subspace 7idctccted formed by the states in which the particle is detected at some time at X , and (1 — P) is the 
projector on the subspace 7i n ever detected of states in which the particle is never detected at X . The fact that those 
two classes of states form orthogonal linear subspaces follows from the superposition principle again. 

Thus, T is properly defined by (||) on 7idctectcd only. If we define the time-of-arrival operator by (^|) on the entire 
state space, then we have the awkward consequence that the states in the range of P(T = 0) and the ones in 
Wnever detected are both annihilated by T. Namely T does not distinguish the states in which the particle is detected 
at T = from the ones in which it is never detected. 

The full information that we need in order to compute 7r(T) is contained in the incomplete spectral family P(T), 
or, equivalently, in the two mutually commuting operators P, and T, where T is a self-adjoint operator on the Hilbert 
space "^detected = Range(P(x)). 

Notice that 

- (iP\Pf(x)P\^) 

is the expected time-of-arrival in those states that are detected at all, and is thus a conditional expectation value. 

As defined in (^), T annihilates all states in 7i aevel detected- It is useful to replace this definition by fixing the 
following convention for the action of T on Hnever detected- We define T on the entire state space, by mimicking what 
happens in classical mechanics: we choose (arbitrarily, at this stage) a (diagonalizable) action of T on TL nevev detected 
with a complex (non-real) spectrum, with the understanding that any complex eigenvalue be interpreted as "the 
particle is never detected". If we use this convention, the operator T is not self-adjoint, but it still has a complete 
and orthogonal basis of eigenstates. The reason for such a convention (which we postulate from now on) will be given 



below; its utility will be particularly clear in |10|. 



C. Partial characterization of T from its classical limit 

In the previous section we have argued on general grounds that a time-of-arrival operator T giving the time-of-arrival 
probability distribution should exist. How can we construct the operator T(x), from the knowledge of the dynamics 
of the system? Let us work in the Heisenberg picture. The quantum theory is defined by the Heisenberg state 
space, in which states do not evolve. Let ip be a Heisenberg state. The elementary operators are Heisenberg position 
operator xq and the momentum operator po, representing position and momentum at t = 0. Since all operators can 
be constructed in terms of xq and po, we expect to be able to express T as an operator function of xq and pQ. A key 
requirement on T is that it yield the correct results in the classical limit (Bohr's correspondence principle). If so, the 
dependence of T on xq and po should reduce to the classical dependence of T on xq and po in the classical limit. This 
indicates that the dependence of the operator T on xq and po is given by some ordering of the function ([!]), which, 
we recall, was obtained by inverting the solutions x(t;xo,Po) of the classical equations of motion. Thus, we should 
have that 

f = t(X;x a ,p a ) (10) 

where an ordering has to be chosen. The c-number X, we recall, is the position of the detector. Notice at this point 
the usefulness of the convention that complex eigenvalues represent non-detection: this can go through naturally in 
the classical limit. Eq. ( |To| ) does not suffice in general for characterizing T uniquely, because the correct physical 
ordering of the operator function can be highly non-trivial. In the companion paper [|10j, we investigate a technique 
for constructing the operator T{X) and fixing the ordering ambiguities. 

In order to emphasize the dependence of the time-of-arrival operator T on the position X of the detector, we will, 
from now on, write the operator as T{X). Analogously, we will denote the spectral family of projectors associated to 
T(X) as P(T;X), and the probability distribution of the time-of-arrival at X as ir(T;X). 

The construction above can be easily generalized to a systems with n degrees of freedom. A classical state of 
such a system is described by a point in the 2n~ dimensional phase space T, with coordinates {z l , i = 1 • • • 2n}. The 
dynamics is generated by the Hamiltonian H(z l ). The Hamilton equations of motion are dz l (t)/dt — {z 4 ,iJ}. The 
general solutions to these equations can be written as 

z\t) = z i {zi,t), i = l---2ra, (11) 
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where Zg are 2n independent integration constants. In particular, we may choose as integration constants z l the 
values of z l at t = 0, z l (0). The above equations enable us to compute the state of the system at any time t. We are 
interested in the time-of-arrival of the particle at a given value of one of the phase space coordinates, say z . (Since 
the coordinates z 1 are arbitrary z 1 can be any combination of dynamical variables.) To compute the time-of-arrival 
T(Z) at z 1 = Z, we solve the 1st equation of the system ( |TT| ) 

z 1 (z \,t) = Z (12) 

with respect to t, obtaining t(Z; Zq). The time-of-arrival T(Z) is then given by 

T(Z) = t(Z;zi). (13) 

Now, in the quantum theory the constants of motion Zq correspond to Heisenberg operators Zq. Eq. ( |l0| ) is 
immediately generalized by "quantizing" (|D|) as 

f(Z) = t(Z;zi). (14) 

where, again, the operator T(Z) is given only up to the ordering. 

Before concluding this section, we add an important comment on the seemingly puzzling case T < 0, namely when 
the detection time is earlier than t = 0. In the classical case, the particle can be detected without being disturbed, 
but not so in quantum mechanics; therefore one might wonder about the meaning of a detection at T < for a particle 
that has a certain state at time t = 0. The difficulty is avoided by choosing the definition of "state" appropriate to 
the present context. Consider the classical case first, xq and pa fix a unique solution of the equations of motion. 
This solution could be characterized by the values of x and p at any other time, or any two constants of the motion. 
The definition of "time-of-arrival" that avoids the problem of detection-before-preparation is the following. We are 
interested in the arrival time of a particle which is moving according to the (unique) solution of the equations of motion 
characterized by the fact that at t = the particle is at xq and po if not disturbed. Analogously in quantum mechanics 
T is the time-of-arrival of a particle that at an earlier time t (arbitrarily in the past) was in the (Schrodinger) state 
ip(t) uniquely characterized by the fact that, if not disturbed, it would evolve to the state "0(0) at i = 0. In the last 
section we shall describe a general way of dealing with this situation. 

To summarize: In this section we have put forward two physical hypotheses: 
• The probability for the time-of-arrival 7r(T; X) -an experimentally measurable quantity- can be computed by 

tTtP(T;X) = (<P\P(T;X)\tP), (15) 

where P(T;X) are the projectors on the real component of the spectrum of a diagonalizable operator T(X). 
The states in the span of the non-real component of the spectrum of T(X) are never detected at X. 



The operator T(X) is given by a suitable choice of ordering from the equation 

f(X)=t(X;x ,p ) (16) 



or, in general, Equation (14) 



The first hypothesis is motivated by our confidence in the general validity of the superposition principle. The second 
hypothesis is motivated by our confidence in the correspondence principle. In the next section we investigate some of 
the implications of these hypotheses and we illustrate the construction and the use of the operator T(X) in a simple 
case. More interesting models, with a non trivial P operator, will be presented in the companion paper p0| . 

III. TIME-OF-ARRIVAL OF A FREE PARTICLE 

Consider a nonrelativistic free particle in one dimension. Dynamics is generated by the Hamiltonian H := p 2 /2m. 
The solutions of the classical equations of motion are 

x(t; x ,Po) = — t + x (17) 

TO 



5 



The inversion of these yields the time at which a particle that at t = has initial position and momentum Xo,f>o * s 
detected at the position X (as in equation ([j])) 

T(X)=t(X; XQ ,po)= m{X ~ Xo) . (18) 

Po 

Notice that up to problems at the point po = (problems with which we shall deal extensively later) the particle is 
always detected. In particular, T(X) is never complex. This simplifies the setting greatly, since we may disregard the 
complications arising from the existence of (finite) regions of phase space in which the particle is not detected. 

Let us consider the usual quantum theory of a free particle. We work in the Heisenberg picture. We have Heisenberg 
(non-evolving) states tp, and time dependent Heisenberg position and momentum operators x(t),p(t), expressed in 
terms of £q andp - Following the ideas of the previous section, we explore the hypothesis that the quantum probability 
distribution 7r(T; X) of the time-of-arrival at X of the particle can be computed in terms of an operator T(X) defined 
by a suitable ordering of the (formal) operator function 

t(x) = "(*-*>) . (19) 

Po 

Notice that the Heisenberg position operator is 

x(t) := e iAt / H x e- i6t / n = ^ t + x ; (20) 

TO 

to be compared with ( p"7| ) : Thought rarely emphasized, classical and quantum dynamics are generically related by the 
equation 

x(t) = x(t;x ,p ) (21) 

where the RHS is an operator function corresponding to an ordering of the solution x(t; xo,po) of the classical equations 
of motion. In general, equation ( pl| ) is of scarce use for solving the quantum dynamics, since the associated ordering 
problem is serious; but in simple cases such as the free particle, we see from (|20| ) that the natural ordering suffices. 

We explore here the possibility that in the case of a free particle a natural ordering suffices for the time-of-arrival 
operator as well. Namely, we study the choice of a symmetric ordering for the operator (|l9|). We thus define, 
tentatively, 

f(X):=^-mJLx JL (22) 
Po VPo VPo 

In order to study this operator, let us choose a concrete representation for the Hilbert space, namely a basis. It is 
convenient to use the momentum basis that diagonalizes po, because this basis makes the definition of — \= simpler. 



Thus, we work in the Heiscnberg-picture momentum basis. The states are represented as functions tp(k) € L 2 (H), 
and the elementary operators Xq and po are given by 

xq o tp(k) = i — ijj(k) and p Q o ip{k) — hk ijj(k). (23) 
dk 

In terms of the above operators, we have, for example, the Heisenberg position operator (|20| ) 

„ / \ Txk d ,^ , . 

i (t) = — t + i—. (24) 
to dk 

In this representation, the operator T(X) given in ( p2| ) is 

TO 1 d 1 TO X 



T(X)o^(k) = 



1 h y/kdk yfk ft k 



Hk)- (25) 



(We always take the principal value of the square root: \f~k~ = i\/\k\ for k < 0.) 

Notice that the 1— parameter family of operators T(X) can be generated unitarily via translations 

f{X) = e- ikx f(0)e ikx . (26) 



G 



Therefore it is sufficient to study the operator T(0), namely we do not loose generality by assuming the detector to 
be at the origin. We thus put X = from now on, and drop the explicit X dependence 



A At N 771 1 <f 1 

T:=T(Q) = -i — -= — -=. 

We will be interested in the operators corresponding to other positions of the detector later on. 
In the momentum representation, the eigenvalue equation for T 



T\T) =T\T) , 



becomes 

fog T (k) 

where we have introduced the notation 



1 d 1 



g T (k) = T g T (k). 



9r{k) 



'k\T) 



(27) 



(28) 



(29) 



(30) 



for the momentum representation of the eigenstate of T. The eigenvalue equation is easily solved (in each half of the 
real line (k ^ 0)) by 



g T (k) = (a+9(k)+a-9(-k)) 



27T771 



Vk 



exp 



2m 



(31) 



where 9(k) is the characteristic function of the positive half of the real line, and a± are constants independent of k. 
In order to fix the relation between a + and a_, let us act on gr(k) by T. A simple calculation shows 



f o g T (k) = T ■ g T (k) - iJ — — —= (a + + iaJ). 



2nh yjfcf 

Thus, in order to satisfy the eigenvalue equation, it is necessary thatf] 



a_ = ia + . 



(32) 



(33) 



At this point, we encounter a difficulty. The operator we have constructed does not have a basis of orthogonal 
eigenstates. This pathology destroys the possibility of giving T the interpretation we want. In the next subsection 
we show that the eigenstates of T are not orthogonal and we discuss a way out from this difficulty. 



A. Difficulties with T and a regulation 

A simple calculation shows that for any two eigenstates of T with eigenvalues T and T' 

^- \ dk* e ^ 2 M = S(T T>) 
2irm J 7r(J — 1 ' 



'TIT') 



dk g T (k)g T '(k) = 



(34) 



The eigenstates fail to be orthogonal. One can also see that f as defined above has no self-adjoint extensions by 
noticing that its deficiency indices are unequal. 

This difficulty stalled us for sometime, and various attempts to circumvent the problem failed. A way out was 
then suggested by Marolf [Q. The idea is to seek an operator that in the classical limit would not reproduce the 
time-of-arrival exactly, but would rather reproduce a quantity arbitrary close to the time-of-arrival. Namely, we want 
to approximate the time-of-arrival with a different quantity, free from pathologies. It is easy to trace the above 



1 Another approach to obtaining this result is to integrate the eigenvalue equation in a small region around k = 0. One then 
obtains the same continuity condition (p3|) on gr(k). 
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pathology to the singular behavior of 1/k at k = 0. Even classically, a state with k = is physically disturbing: either 
the particle is never detected or the particle may stably sit over the detector. Therefore, we seek a small modification 
of ( |l8| ) such that no divergences occur at po = 0. The modified time-of-arrival can perhaps be interpreted as the 
outcome of a measurement by an apparatus arbitrarily similar to a perfect detector, but which does not allow the 
particle to stand still. 

Let us introduce an arbitrary small positive number e. Consider a 1-parameter family of real bounded continuous 
odd functions f t (k) which approach 1/k pointwise. More precisely, we require 



Uk) = - k 



for \k\ > e 
and nonzero for all k =/= 0. 



(35) 



For instance, we may choose 



/«(*) = ~ k 
fe(k) = e- 



for 
for 



\k\>e, 
\k\<e. 



(36) 



Using this, we define the regulated time-of-arrival operator as 



(37) 



to be compared with the unregulated operator (|25|). Notice that on any state with support on |fc| > e the operators 
T e and T are equal. Their action differs only on the component of a state with arbitrary low momentum. As we shall 
see, the probability distribution for the time-of-arrival 7r(T) computed from T e will turn out to be independent of e 
for states with support away from k — -reinforcing the credibility of the regulation procedure we are using. 

Let us study the operator T e . A key point is that T e commutes with Q(k) = sgn(k) — Thus, we can choose a 

basis of solutions of the eigenvalue equation for T e formed by functions of k which have support on positive or negative 
k only. Now T e is a linear differential operator and since f e {k) — > as k — > 0, there is no continuity condition on its 
eigenstates at k = 0. These two related properties lead to a degeneracy in the spectrum. For each eigenvalue T, there 
are two eigenstates, which we choose as having support in the k > 0, k < regions respectively. Namely 



where 



T e \T, 
f e \T, 

[k\T,- 
[ k I T, - 



= r|r,+) 6J 

= T|T,-) eJ 

for k < 0, 
for k > 0. 



We introduce the notation 

e g$(k) = (k\T,±) e 

for the momentum representation of the eigenstates of T e . A simple calculation shows that these are 

— / ~ 1 ( JiT [ k d k\uk'))A. 



t g±{k)=9{±k) 
Explicitly, for \k\ > e we have 



7f=exp i — 

Zirm y/f e \ to 



± c 



^(k) = e(±k) X /^Vke Xp (^-e 2 : 



ihT . 

2m 



(38) 

(39) 
(40) 

(41) 
(42) 



In order to derive various properties of these eigenstates, it is convenient to introduce new coordinates on the right 
and left halves of the real line, as follows. 



zHk) =L dk 'ih 



(43) 
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In the region |fe| > e, Z ± (k) = (k 2 — e 2 )/2. In what follows, we do not need the specific form of Z ± in the region 
|fc| < 0. Note first that in each half, the Jacobian of the coordinate transformation is non-vanishing VTc ^ 0, and thus, 
the new coordinates are strictly monotonic. At the points \k\ — e, Z ± (k = ±e) = respectively. Furthermore, since 
\fe(k)\ — > rapidly enough as |fc| —> 0, we see that both Z ± 6 (-co, oo). In terms of these coordinates, the eigenstates 
are 



e9 ±(Z±) = J J- / .*, e i * T -* k (44) 
V 2-Km ^// e (fe(Z±)) 

and the Hcrmitian inner product between two states ip, <f> is 

/oo 
dZy e (k(Z^(k(Zv))^(k(Z^)), (45) 
'/ -■- "°° 

From (Q,^5|), the orthogonality of the eigenstates is manifest: 

e <T ) 7 ? |T / ,7 ? ') e = ^5(T-T / ). (46) 

In the new coordinates, completeness too is manifest. We can get the same result in the fc-representation with a little 
work: 

£ I"" dT{k\T,r,) ee (T,r,\k') = ^-Y: C dT WW, , ' e'^W-*^)) 



V VfMVfeW) ^ \dZ^/dk\^Uk)^W) 

= 5{k~k'). (47) 
Since it has a complete orthogonal basis of (generalized) eigenstates with real eigenvalues, T e is self-adjoint. 



B. Time-of-arrival probability density 

Following the general theory of section 2, if the particle is in the Heisenberg state ^(k), the probability density w(T) 
of the time-of-arrival is the modulus square of the projection of the state on the T-eigenstates of the time-of-arrival 
operator. Since these are doubly degenerate, we have in the present case 



7r (T) = |e<r, + |^)| 2 + | e <T,-^)| 2 



(48) 



If we assume that the support of tp(k) does not contain (an arbitrary small finite region |fc| < S around) the origin, 
we can choose e < S and, using the explicit form ([42]) of the eigenstates, we obtain the following expression for 7r(T) 



tt(T) 



27rm 



dkVke 



exp 



iTh(k 2 - ( 2 
2m 



+ 



dk\Pk 



exp 



iTTijk 2 - e 2 ) 
2m 



r/>(k) 



(49) 



Notice that the e dependence gives only a phase that disappears when we take the absolute value squared. Namely 



AT) 



2-Km 



dk\fk 



exp 



iThk' 
2m 



dkVk 



exp 



iThk 2 
2m 



(50) 



We thus have the result that for the states that do not include an amplitude for zero velocity, the time-of-arrival 
probability distribution computed (with e sufficiently small) with the regulated operator T c (x) is independent from e. 

The two terms in ( |50| ) correspond to the left and right moving component of the state. Therefore, we have 
immediately that the probability ir + (T) (and tt~(T)) that the particle is detected in X = while moving in the 
positive (or negative) direction is 



TT^T) 



2nm 



dkVk 



exp 



iThk 2 
2m 



i>(k) 



(51) 
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Finally, the result generalizes immediately to the case in which the detector is not placed in the origin, but rather 
in an arbitrary position X . The eigenstates of the operator T e (X) 



f e (X)\T,±;X) e = T\T,±;X) e (52) 
are obtained using the unitary translation operator 

\T,±;X) e = e-^ x \T,±) e . (53) 

yielding 

e9% x (k) = < k | T, ±; X) e = e- lkx ^(k). (54) 

The projectors considered in section 2 are given by 

P(T;X) = \T,+;X) ee (T,+;X\ + \T,-;X) £e (T,-,X\. (55) 

The probability density of being detected at time T by a detector in X that detects particles traveling with positive 
(negative) velocity is 



ft f ±0 ° ( iTTik 2 \ 1 

• dkVkexp I ikX I ibik) 

o V 2m J 



27TTO 



(56) 



Equation ( p6| ) represents our final result for the probability distribution of the time-of-arrival at X of a free quantum 
particle. 

C. Time-of-arrival of a Gaussian wave packet 

As an example of an application of the above result, we compute the probability distribution for the time-of-arrival 
of a Gaussian wave packet. Consider a Gaussian wavepacket localized about a point (say) to the left of the origin at 
time t = 0, and moving (say) to the right. In the standard Schrodinger-picture position representation, let this wave 
packet be given by the following normalized solution of the Schrodinger equation 

Mx t) ( 5 -\ " ^ exp f (^o+^- *o)) 2 ) (57) 



2m 



Expectation values are as follows 



(p(t))=hk and (Ap(t)) = ^ (58) 



t 2 h 2 

(x(t)) = x + Ukot/m and (Ax(t)) =S\ 1 + 9 . (59) 

y 4d 4 m z 

If we choose \xq\ >> 6, \ko\S >> 1, xq < and fco > 0, this state represents a particle well localized to the left of the 
origin and with a well-defined positive momentum at time t = 0. In the Hcisenberg-picture momentum representation 
(P3|), this state is given by 

^{k)=[ — \ cxp(-(k- k Q ) 2 6 2 -ikx ) . (60) 

The envelope of this wave function is a Gaussian of width 1 /S centered at fco . 

Using the theory developed, we can compute the projection of this state on the eigenstates of the time-of-arrival 
operator. We assume here that e can be taken arbitrarily close to 0. (See jl3| for the relevant integrals). We obtain 

ft / 5 2 n 1/4 



5 <T,+;Jf|^> = \/-(^3 ) e- d fc «r(5/4)x 
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exp 



iKT 
2m 
f(2S 2 k 



-5/4 



(26 2 k + i(X -x Q )) ■* 
z(A-xo)) 2 



V 4<5 2 + 2ihT/m 



C 



1 

2' 



5 3 (2S 2 k + i(X -x )) 2 
4' 2' 4(5 2 + 2ihT/m 
(2<5 2 fc + «(^-xo)) 2 



A5 2 + 2ihT/m 



(61) 



IT,- 



\X\i>) = \ - 
m 




x exp 



2 5 7T J 

i/iT 



1/4 



e -* fcs r(5/4)(-*)x 



-5/4 



* +-5— 1 -(2(Tfco + i(X-xo))-$ 



2m 

/ (2J 2 fc + t (X- a:o )) 2 
V 4<5 2 + 2ihT/m 



C 



1 

2' 



5 3 (2J 2 fc +i(X-x )) 2 



4' 2' 4<5 2 + 2ar/m 
(2S 2 k + i(X -x Q )) 2 ~ 
AS 2 + 2ihT/m 



(62) 



where we have ignored the components in |fc| < e, since e can be taken arbitrarily small. m;z] — iFi[n;m; z] 

is the Kummcr confluent hypergcometric function, C[n,a;z] — C^(z) is the nth generalized Laguerre polynomial and 
r(n) is the Euler gamma function.) The probability distribution is then given by (^). The expression above for the 
probability distribution of the time of arrival of a Gaussian wave packet is a bit heavy; in order to unravel its content, 
we have expanded it in powers of small quantities in the Appendix A, and we have plotted the total probability 
density 7r(T) at various detector positions X in Figure 1 (choosing a Gaussian state (|57| ) with x — —5, fc = 20, 8 = 
.5, h = m = 1). To begin with, the term corresponding to negative velocities is exponentially small; indeed, it derives 
from the scalar product of a Gaussian wave packet concentrated around a positive k with a function having support 
on k < 0. 




FIG. 1. Time-of-arrival probability densities tt(T; X), plotted at the detector positions X — —5, —3, — 1, 1, 3, 5 . 



The total detection probability density at time T for the detector in position X is a function (more precisely, it 
is a density in T) on the T — X plane. This function is concentrated around the classical trajectory of the particle 
X = xq + poT/m, with a (quantum) spread in T that increases with the spread of the wave packet, namely with the 
distance of the detector from the initial state. 

In Appendix A, we compare our result with the probability density obtained indirectly using the Shrodinger proba- 
bility density. We find good agreement within second order in the deBroglie wavelength of the particle, and we discuss 
the order of the discrepancies. Thus, our result is reasonable to leading order. Whether or not it is physically correct 
to all orders is a question that can perhaps be decided experimentally. A discrepancy with an experimental result 
may indicate an incorrect ordering of the time-of-arrival operator, or a more general difficulty with our approach. 



11 



IV. DISCUSSION 



A. Time representation 



Anytime we have a self-adjoint operator in quantum mechanics, we may define a representation that diagonalises 
this operator. Namely, we may choose the eigenbasis of the operator as our working basis on the theory's Hilbert space. 
Nothing prevents us from doing so with the operators T e (X) as well. Let us therefore introduce a "time-of-arrival 
representation", or, for short, a "time-representation". We fix an e and define 



1 pt(T > X):= e (T,±;X\i>)= J dk e g% x (k) ^(k). (63) 

Clearly, we can do quantum mechanics in the (T, X) representation, as well as we do quantum mechanics in the 
position, momentum, or energy representations. Since the eigenstates e g^ x (k) have support on positive k, we must 
interpret ipj(T,X) as the amplitude for the particle to be detected by a detector placed at X in an infinitesimal 
neighborhood of T coming from the left, and ^(T, X) as the amplitude to be detected at A in a neighborhood of T 
coming from the right. 

What is the relation between the ip^(T,X) amplitude and the conventional Schrodinger wave function 2p(x,t)l 
Notice that the first is defined by ipf(T,X) = e (T,±;X\tp), where |T,±;A) e is an eigenstate of T e (X) with 
eigenvalue T; while the second can be viewed as defined by ijj{x,t) = (x;t\ip), where \x;t) is the eigenstate of the 
operator x(t) with eigenvalue x. At first sight, the two seem to be related to the same quantity (up to the e and 
the distinction between the two directions of the velocity): they both refer to probabilities of being detected at some 
space-point and at some time. However, this naive observation is very misleading. The quantity \i/>(x, t)\ 2 dx is the 
probability in space that the particle happens to be between the positions x and x + dx at time t, as opposed to being 
elsewhere at time t. Whereas, the quantity \tp (T, X)\ 2 dT is the probability in time that the particle happens to 
arrive between times T and T + dT at the position A, as opposed to reaching X at some other time. The two bases 
| T, ±; A ) e and \x;t) are two well defined (generalized one-parameter families of) bases in the Hilbert space, but they 
are distinct. 

In particular, the two bases | T, ±; A ) e and | x; t ) have distinct dimensions, because \ip(x, t)\ 2 dx and |V' ± (T, A)| 2 dT 
must both be dimensionless probabilities. Thus, the transformation factor between \T, ±;A) and \x;t) has the 
dimension of the square root of a velocity. Indeed, let us write the two (generalized) states explicitly in the Heisenberg 
momentum representation. Restricting ourselves to k > and taking e to zero for simplicity, we have from ([42] ) 

(k\T,+;X) e = J±JI eX p ( % ^k 2 - t kx) , (64) 
V 27r V m \ 2m / 



while, as it is well known, 



[k\x,t) = J-^-exp( ^k 2 -ikx) . (65) 



2tt V 2m 



Therefore, taking x — X and t — T we have 



{k\T,+;X) £ = J— {k\x,t). (66) 
V m 



A physical understanding of the curious J ^ factor that characterizes the eigenstates of the time-of-arrival operator 

can be obtained as follows. Consider a well localized wave packet travelling with velocity v = ko/m. We have 
approximately |^ ± (T, A)| 2 w v \ili(x,t)\ 2 . Now, consider the x — t plane. The wave function ip(x,t) is significantly 
different from zero on a band around the classical trajectory. The classical trajectory is a straight line with a slope 
given by the velocity v. The ratio between a vertical and an horizontal section of the band is therefore precisely 
v. Thus, in order to have both total probabilities normalized to 1 when integrating along at = constant, or a x — 
constant line, the probability density in space and the probability density in time must be related by a factor v. 

In Figure 2 we have plotted the usual Schrodinger probability density in x - ir(x;t) = \t/j(x,t)\ 2 - at various times, 
for the same state used for Figure 1. 
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FIG. 2. Schrodinger probability densities n(x;t) for position, plotted at the times t = 0, .1, .2, .3, .4, .5 



The time-of- arrival probability densities are v — ko/m times as large as the Schrodinger densities for values of 
X = x and T = t near the classical trajectory. 



B. Time-energy uncertainty relation 

Position-momentum uncertainty relations are of wide use in quantum mechanics, and can be cleanly derived from 
the formalism on very general grounds. The commutation relations [x,p] — —ihl imply AxAp > 2h. Time-energy 
uncertainty relations are of wide use as well (e.g. between the width of a spectral line and the lifetime), but their 
general derivation is notoriously tricky. If one tries to reproduce the position-momentum derivation for the time-energy 
case by assuming the existence of an operator T such that 

[f , H] = -ihl (67) 

from which 

ATAE > 2h (68) 

would follow, then one would clash against a well known non-existence theorem for T. The theorem states that the 
commutation relations [A, B] — —ihl between two self-adjoint operators A and B implies that the spectrum of both 
operators is the real line. However, the spectrum of the Hamiltonian is bounded from below in all reasonable systems. 
Ergo a time operator T satisfying (p?]) does not exist. 

This theorem might have been the reason for which a time-of-arrival operator has virtually never been considered in 
quantum mechanics. In fact, it is often stated that time cannot be an operator in quantum mechanics, with the above 
theorem as a proof. Here, we show how one can rigorously derive time-energy uncertainty relations for a quantum 
particle, and how the existence of T e circumvents the theorem. 

The commutation relations between T e and the Hamiltonian are easy to compute. In the momentum representation 
we have 

[f e ,H] = -ihi-ihh e (k), (69) 

where 

h e (k) = l-kf e {k). (70) 

The function h e (k) is bounded (by e 2 , if we chose f e (k) as in (|36|), which we do here for simplicity) and has support 
on the small interval \k\ < e. For the particle in the state ip(k) the following uncertainty relations follow 



AT e AE > 2h — Aire 2 J dk \ip{k)\ 2 > 2h - Aire 2 . (71) 
By chosing e sufficiently small, we obtain an uncertainty relation that approaches (pq) to any desired precision. 



13 



C. On the definition of state in quantum mechanics 



Finally, let us return to the problem we briefly discussed at the end of section 2, which is the interpretation of the 
time-of-arrival T when T < 0, namely when detection time is earlier than the time t — at which the initial state is 
given. We have suggested that in this case the correct interpretation of T is the following. T is the detection time for 
a state that arbitrarily in the past was in a state that would have evolved to the t = initial state if undisturbed. A 
cleaner way of dealing with the general situation, is to make use of a fully time-independent notion of "state" and a 
fully time-independent version of phase space and quantum state space. This can be done as follows. Consider first 
classical mechanics. Let us denote a single solution of the equation of motion as a "physical history" of the system. 
(A physical history should not be confused with the histories considered in sum-over-histories theories: a "physical 
history" here is a history satisfying the equations of motion.) Let V s be the space of these physical histories. A point 
in T s represents an entire evolution of the system. r s can be coordinatized by the 2n integration constants z % . We 
ask for the time-of-arrival at a; of a system following one of the motions in Y s . This time-of-arrival is given by (|]1]) . 
The key to the matter is that there is no need to choose a time in order to specify a physical history. 

In quantum mechanics we can define the Hilbert space TLs of the solutions, of the Schrodinger equation. A vector in 
Hg represents an entire (quantum) motion of the system, without reference to any particular time. The conventional 
Heisenberg operators are defined on Hs- The operator (|lj) is properly defined on Hs- We may choose to represent 
the vectors in Hs by means of the value that the Schrodinger state would take (if undisturbed) at t = 0; therefore it 
makes sense to define T as a function of the operators xo and po which are defined on the states at t = 0. 

In this regard, it is interesting to notice that the original definition of the "Heisenberg picture Hilbert space" given 
by Dirac in the first edition of " Principles of Quantum Mechanics" is the definition of Hs given above fTl j. It is only 
later that the "Heisenberg picture Hilbert space" came to be mostly identified with the state space at a fixed time 
(both interpretations of Tig can be found in the literature). A crucial advantage of using the definition of Hs given 
here is that this definition can be extended to systems without Newtonian time at all |l4|,[l6|]. We will exploit this 
point of view in jlG] . 



V. CONCLUSIONS: X «-» T EQUIVALENCE IN QUANTUM THEORY 

Let us summarize our results. We have considered the problem of computing the time-of-arrival T of a quantum 
particle at a position X. Relying on the general validity of the superposition principle, we have argued that the 
probability distribution for T can be obtained by means of an operator T. This operator is, in general, not self- 
adjoint. However, it admits an orthogonal basis of eigenstates. The eigenstates with real eigenvalues correspond to 
(generalized) states for which the detection time is sharp. The eigenstates with complex eigenvalues correspond to 
states that are never detected. 

The time-of-arrival operator is partially characterized by its classical limit, which fixes its dependence on the 
position and momenta operators, up to ordering. We have considered the simple case of a nonrelativistic free particle, 
using a tentative natural ordering. A regulation procedure allows us to to find a self-adjoint time-of-arrival operator 
(there is no non-detection in this case), and we have studied the probabilities the operator yields. We showed that the 
regulated time-of-arrival operator can be used to derive time-energy uncertainty relations, circumventing a well-known 
nonexistence theorem, and that it yields a well-defined "time representation" for the system. 

For more general systems, the classical limit is not likely to be sufficient for constructing the operator. In a 
forthcoming companion paper [^p| we investigate a general technique for constructing the time-of-arrival operator in 
general cases. We will study in detail the particle in an exponential potential, where the P operator is non-trivial. We 
will also investigate parametrized systems and theories without a Newtonian time, arguing that the ideas presented 
here may be relevant for the interpretations of quantum gravity. 

We close with a general comment on the significance of the result obtained. In classical mechanics there is a 
hidden equivalence between the independent time variable t and the dependent dynamical variables -the position x 
in the present case. This equivalence is made manifest by expressing the theory in parametrized form, namely by 
representing the evolution in terms of the functions t(r) and x(t) of an arbitrary parameter r. In more elegant and 
general mathematical terms, there are formulations of mechanics, e.g. the presymplectic formulation, in which the 
distinction between dependent and independent variables is inessential; see for instance Arnold's classic text [ p"7[ . 
A parametrized representation is commonplace in special relativity (where T is called x ) since it allows manifest 
Lorentz covariance. 

It is commonly stated that this equivalence between independent (time) and dependent variables is lost is quantum 
mechanics. The arguments in support of this claim are common in the quantum gravity literature and take various 
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forms. For instance, it is said that the wave function must be normalized by integrating in x and cannot be normalized 
by integrating in t. Or it is said that probabilities are always probabilities of different outcomes happening at the 
same time, never at the same position. It is our impression that these claims are misleading. The mistake is to assume 
that the x <-> t equivalence has to be realized as an equivalence in the arguments of the Schrodinger wave function 
ip(x,t). 

The conventional formulation of quantum mechanics in terms of the Schrodinger wave function ip(x,t) has already 
broken the x <-» t equivalence. Indeed, it is a formulation tailored to answer the (experimental) question: "What 
is the probability of the particle being here now, as opposed to that of being elsewhere now?". The corresponding 
representation diagonalises the Heisenberg operators x(t). It is the experimental question considered, and the related 
choice of basis, that breaks the x <-> t equivalence. Quantum mechanics allows us to consider the following question 
as well: "What is the probability of the particle getting here now as opposed to getting here elsewhen?" . In order to 
answer this question, one is led naturally to the ^ ± (T, X) representation in which the role of position and time are 
to a large extent interchanged. In particular, the wave function is normalized in time and probabilities of events at 
the same position are considered. 

To avoid misunderstandings, let us make clear that we certainly do not claim that space and time have the same 
nature, nor that their role in the quantum mechanics of a particle is exactly the same. What we suggest is that the 
common arguments that x and t can be treated on equal footing in classical mechanics, but not in quantum mechanics, 
might loose force under closer scrutiny. Contrary to the above arguments, our analysis has revealed an underlying 
hidden equivalence between "dependent" and "independent" variables in the quantum theory of a free particle. 
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APPENDIX A: IS THE COMPUTED PROBABILITY DENSITY REASONABLE ? 



Let us now investigate whether the result we have obtained for the Newtonian free particle is physically reasonable. 
The simplest check is to compute the expectation value of the time-of-arrival of a wave packet with initial position 
xq and momentum po. The result should agree with the expected classical time-of-arrival T = t(X;xo,po)- Since 
the operator was constructed via a factor-ordering and regulation of the classical solution, the expectation values 
obviously satisfy Ehrenfest's theorem. A more accurate check which goes beyond the semiclassical approximation 
is to compare the probability distribution we have obtained with the one we can estimate by indirect but intuitive 
methods. In appendix A.l, we first compute an approximation to the probability amplitude ( |6l| ) -and the resulting 
probability density- for a Gaussian state. This also gives us some intuition for the behavior of the the distribution. 
Then in appendix A. 2, we compute the Schrodinger current through the detector position, and compare this with the 
approximation we will obtain in section A.l. 



1. Gaussian state 



Consider the Gaussian wavepacket of section 3.5. In the momentum representation, this state is given by 



2S 



2 \ 4 



ip(k) = I exp (— (k - ko) 2 8 2 - ikxo) 



(A.l) 



The envelope of this wave function is a Gaussian of width 1/8 at positive momentum fco. In order to simplify the 
calculations, we slightly modify this state by assuming it to be zero for k < e. Clearly, since this modification is far 
out on the tail of the Gaussian, the error we make is very small (More precisely, one can show that it vanishes as 
exp (— 1/koS)). We thus replace (AT) by 



25 2 \ 3 

ip(k) =6{k-e)[ exp (-(& - k Q ) 2 5 2 - ikx ) 



(A.2) 



where 9{x) is the characteristic function of the positive line. Substituting this expression into (|6^), we find that the 
amplitude for the particle to be detected between T and T + dT at the position x is 



i>+(T,X) 



2 Ti m \ 7T 



OX2 \ POO _ / +rp \ 

—J J dkVkexp\-{k-k ) 2 S 2 ~i — k 2 +ik(X -x )\ 



2S 2 



2\ 1/4 r 



dkVkexp (-(k - k ) 2 S 2 - i—k 2 + ik(X - x )j - i?i 



(A.3) 



where R\ is again exponentially small, ~ exp(— l/kgS), and we have disregarded an e-dependent phase, which is not 
going to affect the probability distribution. 

Now we can expand \fk about the peak of the Gaussian as 



Vk- 



p- I k~k Q . n fk-k 



(A.4) 



Using this expansion, we can compute the Gaussian integral in (A.3) 



dkVkcxp ( — (k — ko) 2 S 
\ 



2m 



_ ik (X-x a )-i^kl 



jfegTT 



2m 



exp 



(X-xo-^) 2 
^ 2 + iE) 



l + i 



X - Xq 



hTkg 



u v 2m 



(A.5) 



where 



\Ba\ < 



8A-, 



3/2 



dkk e 



2-k 2 6 2 



^QTT 



1 1 



16k 2 S 2 S 



(A.6) 
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is of second order in Thus, to first order in the (small) quantity j^, we have 

* ( ' )_ V mUJ 4(^ + i^)V a ^ 1 + '4* ft (^ + i^)J' 

2 



To this approximation, the probability density 7r+(T;x) = |^ + (T;x)| 2 is then 

Due to the Gaussian factor, the probability distribution is centered on the classical time-of-arrival T — m X f ~ k ^ a , with 
width 5, and vanishes exponentially outside such a region. 

Notice also that due to the Gaussian factor the third term in the numerator of the amplitude of the Gaussian 
is of order (j^) 2 in all the region where the probability density is not exponentially small, so we may rewrite the 
probability density as 

as an approximation correct to order {jj^)- 

2. Current 

In ordinary quantum mechanics, one can define a current whose time and space components are 

j°(x,t)=p(x,t)=il)(x,t)^(x,t) (A.10) 

j\x,t) = (A. 11) 
zmi 

where il>(x, t) is the quantum state in the conventional Schrodinger-picture position representation. Since the state 
satisfies Schrodinger's equation, this current is conserved 

dof + dd^O. (A. 12) 

Consider, for a particle in one dimension, the spacetime region (a half strip) defined by x < X, t\ < t < t2, and 
integrate the divergence free current over this volume. Dropping the boundary terms at x = — oo, we find that the 
outgoing (i.e., rightward) flux of this current through the timelike boundary at x = X in the time interval (ii,*2) is 



rti 

Jt! -"" ' 



ft f — — I 

dt -^ = ^~ dt (ipd x ip - ipd x ip)\ x=x 

/X 1 r-X 

dxp\ tl - / dxp\ t2 . (A.13) 
-oo J — oo 

The last line in the above equation is the probability that the particle is found in the LHS region (X < x) at time t\ 
minus the probability that the particle is found in the LHS region at a later time ti. It is tempting to identify the 
flux between the times T and T + dT (through the timclikc surface at x — X) of the current density 

jx(T)dT = Aj m (jfaT)d x il>(x,T)) dT (A.14) 

as the probability density that a detector placed at x = X will detect the particle. Note that the Schrodinger current 
density has the correct dimensions of a density in time, namely [T] -1 . The problem, of course, is that the current 
represents the net flux of probability across x = X, and thus corresponds to the difference between the probability of 
crossing to the right and the probability of crossing to the left. Namely we may expect that 
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j x (T)~ir+(T;X)-n-(T;X). 



(A.15) 



within some approximation. In fact, the current jx(T) is not positive definite, and we believe that this is related 
to the difficulties in the approaches of Q and [||. Equivalently, we may try to interpret the current of a pure right 
moving state as the probability density that the particle crosses x = X. We are reassured in doing this by the fact 
that in this case the current is positive definite, and its integral over all times gives one. This conclusion is based on 
the assumption that the particle cannot "zigzag" across the X = x line, an assumption which might be valid only if 
we look at sufficiently large times. 

So, when ip is a pure right moving state, the rightward flux density is positive^] and integrates to 1 over all time. 
(To see this, take the limits t\ — > — oo and t-2 — > oo in (A. 13).) It is therefore at least consistent to interpret this 
as an estimate of the probability density. Does this estimate yield the correct semi-classical limit? The expectation 
values of the usual position and momentum operators satisfy Ehrenfest's theorem, since the probability densities are 
associated with the decomposition of a state onto the spectrum of some sclf-adjoint operator. The same is true of the 
probability densities we have obtained from the time-of-arrival operator. The above current density is not obtained 
via a spectral projection, however, and is not associated with a "time operator" . How do the expectation values of 
the time-of-arrival behave wrt. this estimated probability density, and do they correspond to the classical limit in 
some way? We next proceed to analyze this issue. 

The expectation value of the time-of-arrival of the particle at the position X is naturally defined only when the 
state has support only on the k > region and is then given by 
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(A.16) 



where we have dropped a surface term at k — since we assume that ip(k)k^o ~ * k 2+ . The final expression in the 
above equation is the symmetric factor ordering of the classical expression T = m(X — xq)/pq, thus we do recover the 
correct semiclassical limit. 

For the localized right moving wave packet previously considered, the current is easily computed, giving 
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(A.17) 



The integral can be explicitly done, yielding 
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(A.18) 



which is precisely the approximate form for the probability we computed from the time-of-arrival operator (see (A. 9)). 
Thus, the probability distributions computed with the time-of-arrival operator and by means of the current agree for 
a right moving localized wave packet to order (l/k S), which is one order beyond the classical limit. 

For a general state, roughly localized in momentum state around a momentum kg, we can compare the current 
(A. 14) with the difference between the probability of being detected moving towards the right minus the probability 
of being detected moving left. Namely we can estimate 



See (A. 14) and recall that this is a free particle. 
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d x {T)=j x (T)-{-K+{T)-v-{T)). 



(A.19) 



Explicitly, using ( A. 14 ) and (|50|), we have 

dx(T) = —Im U(X;T)dxi>(X;T)) 
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We now put X = for simplicity, without loss of generality. A little algebra gives 
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(A.21) 



Now suppose that our measurement has a precision ST. Then the difference of the probability densities, averaged 
around T is 
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The integral in dT' can be done easily, which yields 
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(A.23) 



For large ST, compared to the "deBroglie time" of the particle |£r (where k is the highest momentum in the support 
of the wave function), the factor { } in the integrand approaches a delta function in the difference between the two 
integration variables (plus a term with a delta function in the sum of the two integration variables, which we may 
assume to be negligible for states sufficiently localized in momentum space), and the integral is then suppressed by 
the (y/k — \fk!) 2 factor. This indicates that the two ways of computing the probability for the time-of-arrival approach 
each other when our resolution time is larger than the particle's deBroglie time. 
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